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ABSTRACT 

We present results from a Hubble Space Telescope ( HST ) program characterizing the atmospheres 
of the outer two planets, in the HR8799 system. The images were taken over 15 orbits in three near- 
infrared medium-band filters - F098M, F127M and F139M - using the Wide Field Camera 3. One 
of the three filters is sensitive to water absorption band inaccessible from ground-based observations, 
providing a unique probe of the thermal emission from the atmospheres of these young giant planets. 

The observations were taken at 30 different spacecraft rolls to enable angular differential imaging, 
and the full data set was analyzed with the Karhunen-Loeve Image Projection (KLIP) routine, an 
advanced image processing algorithm adapted to work with HST data. To achieve the required high 
contrast at sub arcsecond resolution, we utilized the pointing accuracy of HST in combination with 
an improved pipeline designed to combine the dithered, angular differential imaging data with an 
algorithm designed to both improve the image resolution and accurately measure the photometry. 

The results include F127M (J) detections of the outer planets, HR8799 b and c and the first detection 
of HR8799b in the water-band (F139M) filter. The F127M photometry for HR8799c agrees well with 
fitted atmospheric models resolving a long standing difficulty to model the near-IR flux for the planet 
consistently. 

Subject headings: brown dwarfs - planetary systems — stars: atmospheres — stars: low-mass 


1. INTRODUCTION 

Within the past decade, direct imaging of exoplanets 
has advanced from initial discoveries of favorable low con- 


trast, planetary mass companio ns like 2M1207b (Ch au- 
vin et al |2005a ) and AB Pic b ( Chauvin et al.|20C 515) to 


the HR8799 planetary s 1 

astern (Marois et al.|2008 2010), 

0 Pic b (Lagrange et al. 

2010), and HD 95086b (Rameau 

et al. 2013), and GJ 50 

4b (Kuzuhara et al. 

20T3). Di- 

rectly imaged exoplanets torm a critical su 

bset ol the 


exoplanet population for which it is possible to char¬ 
acterize atmospheric properties from thermal emission. 
Given the stee p, mo n otonic declin e in planet brightness 
with time (e.g. Burrows et al.|2001), the currently known 


directly imaged exoplanets (companions to young stars), 
enable investigations into the early evolution of exoplan¬ 
ets. Imaged young atmospheres als o present a val uable 
comparison to the transmission(e.g. IDeming_et_^LJ[2005 


Sing et aT1|2008|) and emi ssion (e.g. [Charbonneau et al. 
2005||Knutson et al.|2007 ) spectra of intensely irradiated 


planets orbiting older stars. 

HST has been employed successfully in characteriz¬ 
ing exoplanet systems for over a d ecade ([Charbonneau 
et al. 2002 Soummer et al. 2012). We present Wide 
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Field Camera 3 imaging of the multiple planet system, 
HR8799. The system has four you ng, massive plan¬ 


ets (~ 4-7 Mj up ) orbiting a ^30 Myr (Zuckerman et al. 


20111 A5 star, amenable to both photometric and spec- 


troscopic follow-up. This is the first time that planets 
with high contrasts have been successfully imaged with 
WFC3/IR, a camera without a coronagraph and under¬ 
sampled point spread function (PSF). The techniques de¬ 
scribed in the paper should enable future studies of plan¬ 
etary systems with HST. In this paper, we characterize 
the outer two planets, HR8799 b and c, by combining 
the new space-based near-infrared photometric data to 
the existing suite of ground-based data available for the 
two planets and discuss the implications of the new pho¬ 
tometry. 


2. OBSERVATIONS 

Observations for this stud}|^] were taken from Nov 
2011 to Dec 2012 with the near-IR channel of WFC3 
on the HST. Due to the brightness of the primary star 
(J=5.3 mag), all combinations of sub-array and detector 
readout time results in saturated zeroth read. In the ab¬ 
sence of a coronagraph we designed the observations to 
saturate the central ^0."5 of the stellar PSF using 2.2s 
long observations, focusing the program on the b and c 
planets. 

The HR8799 system was observed in three medium- 
band filters, F098M, F127M, and F139M over 15 orbits. 
The orbits were grouped contiguously in sets of three to 
increase PSF correlation. The observing sequence was 
designed to maximize the total roll angle and minimize 
variations in the PSF of the observations over the dura¬ 
tion of the program. To maximize the rotation, the star 

8 GO Program #12511 

























































F098M 


FI27M 


FI39M 


HR 8799b 

■ 

HR 8799b 


HI) 

HR 8799b 




HR 8799c 

1 HR 8799c 1 


HR 8799c 


(fD 




1 ■ 





F098M 

FI27M 


| FI39M 


Fig. 1.— Multi-wavelength HST HR8799 data. Top: Image zones with locally optimized KLIP detections of HR8799b in all three filters 
- F098M, F127M and F139M, respectively. Bottom: Same as above, for HR8799c. HR8799c is only detected in the F127M filter, with 
upper limits in F098M and F139M. 


was observed over two different roll angles in each or¬ 
bit, resulting in observations taken over 30 separate roll 
angles covering ~270 degrees (not uniformly sampled) 
to enable A ngular Differential Imaging (ADI) reduction 
(Marois et al. 2006). Within each roll, the telescope 
was dithered using a customized nine-point spiral dither 
pattern covering a O/'ISxO/'IS region, with half-pixel 
dithers (~64 mas), ensuring a well sampled PSF. This 
observing procedure enhanced stability, but resulted in 
reduced total exposure time on the object and a limited 
range of roll angles within each block of three orbits. At 
the conclusion of the observing program, the full dataset 
comprised 270 images in each of the three filters. 


3. DATA CALIBRATION 

The WFC3/IR images are undersampled data at our 
wavelengths of interest. The DRIZZLEPAC software is the 
standard tool to improve the resolution of HST data, but 
it not optimal for high contrast imaging since it corrects 
the geometric distortion prior to improving the resolution 
of the images, the opposite being more appropr iate h ere. 
We therefore modified the AL ICE pipeline (Soummer 
et al.|2014 Choquet et al. 2014) to include the calibration 
of dithered VVKC3 images. 

The program used a 9-pt spiral dither pattern with 
0.5 pixel dither steps. HST pointing stability within a 
single orbit for small dithers is ~2-5 mas. Using the 
WCS information and the “interleave method” (Lauer 


1999), the dithered frames were combined to improve 


the resolution by a factor of 2 (~64 mas/pixel). The se¬ 
lected dither pattern does not permit reconstruction of 
two fully independent images using interleaving and one 
dithered frame has to be shared between them. The addi- 



Fig. 2.— 5-<7 HR8799 detection sensitivity with AF127M mag¬ 
nitudes as a function of angular separation from the star. The 
black circles are our photometry values for HR8799bc. The solid 
red line shows the contrast estimated from a global KLIP reduc¬ 
tion and the red dashed segments are contrast levels obtained in 
the optimized zones. This figure is intended to help preparation 
of future observations at high-contrast with HST/WFC3 with the 
green line indicating the photon noise floor for the full dataset. 
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tional time-sampling provided by generating two images 
per 9-pt dither meant a deeper contrast and was thus 
preferred over the one image per 9-pt dither case. The 
upsampled images are then corrected from the detector 
distorted frame to an undistorted frame using correction 
maps (J. Anderson, private comm). Finally, the images 
were aligned by cross-correlating on the diffraction spikes 
in the data, resulting in 60 aligned and upsampled im¬ 
ages. 

The data were reduced by the Karhunen-Loeve Image 
Projection (KLIP) algorith m based on prin cipal compo¬ 
nent analysis (PCA, Soummer et al. 2012) used in the 
ALICE pipeline. For each of the 60 upsampled images, 
the reference PSF library is assembled imposing a mini¬ 
mum rotation of 2xFWHM at planet location to reduce 
self-subtraction. 

We assumed knowledge of the location of the planets 
and reduced data in local zones around the expected po- 
sitions for planets b and c, similar to what was done in 
Soummer et al. (2011). The reduction was performed in 
annular sections, over a parameter space exploring two 
radial sizes and three azimuthal widths. With ~25 KLIP 
reductions per image our parameter space includes ~150 
images, giving ^9000 images for each reduction of a given 
planet/filter combination. We explored a number of ge¬ 
ometries and in particular the location of the planet with 
respect to the zone with best results when the planet is 
close to the inner edge of the zone. 

Each reduced image is corrected for the KLIP through¬ 
put loss using forward modelin g, estimate d by projecting 
a TinyTim model PSF |KristjA_al. 20111 onto the KL 
modes, as d escribed in |Soummer et al. (20121; Pueyo 


et al. (2015). We then compute a data quality criterion 


as the ratio of the algorithm throughput to a noise esti¬ 
mate with matched-filtering in a local region around the 
planet. This criterion is therefore proportional to the 
true SNR, but without introducing any of the planet sig¬ 
nal in the calculation so that it does not bias the results 
by amplifying speckles. The reduced data are ordered 
according to this criterion and cumulatively-combined 
using a median. The data quality criterion is then re¬ 
calculated to determine the optimal number of images in 
the final image, shown in Fig. |Tj This approach permits 
exploration of any algorithm parameter space to produce 
a single final image automatically by identifying the best 
combination of algorithm throughput and speckle sup¬ 
pression. 

Photometry was obtained using matched-filtering with 
a truncated model PSF (5 pix) to reduce potential con¬ 
tamination from local speckles. This matched filter com¬ 
bined with partial truncation of the planet PSF by the re¬ 
duction zone leads to an incomplete fraction of the planet 
flux (~ 80%). All these effects were carefully calibrated 
using aperture photometry on archival WFC3 data of 
white dwarfs (GD153, G191B2B) from HST calibration 
programs, and of similar spectral type isolated brown 
dwarfs. The overall photometric correction precision is 
of the order of 3-4%, significantly smaller than the final 
photometric error bars on the planets. 

The PSF library is not rigorously free of companion 
contribution, however the contribution of the companion 
at other roll angles is limited, since the data is not a true 
ADI sequence (the orients span the 270 degrees over mul¬ 
tiple epochs) and the planets are very faint compared to 


the PSF wings. In addition, the noise is zero-mean within 
the reduction zone from KLIP but is not necessarily zero- 
mean within the matched-filter equivalent aperture, and 
small biases may remain. The validity of the forward 
modeling approach for this particular dataset was inves¬ 
tigated by injecting TinyTim PSFs in the data at the 
planets’ radial separation over a range of azimuthal an¬ 
gles (25 PSFs for b and 20 PSFs for c). The injected 
PSF flux was iteratively adjusted to result in the same 
SNR as each true planet after reduction. The mean pho¬ 
tometric error on the synthetic planets was used as a bias 
correction to account for the two effects discussed above 
and the error on the photometry was estimated from the 
standard deviation in the measured signal for each of the 
fakes. Upper limit detections for HR 8799 c in F098M 
and F139M were obtained by adjusting the injected PSFs 
to detect almost all of them. To ensure the absence 
of speckle amplification we injected “zero-flux synthetic 
planets” using the exact same pipeline. No significant de¬ 
tection could be noticed at the location of these injected 
zero-flux PSFs. Various combinations of reduced images 
were tested for the presence of flux at given positions as a 
fun ction of telescop e roll a ngle and reduction parameters, 


Soummer et al. (2011|) for a detailed description. We 


see_ 

estimated the False Positive Fraction (FPF) of 7.6 x 10~ 22 
(SNR ~10) and 3.8xl0~ 7 (SNR~5.3) for b and c in the 
F127M filter assuming Gaussian statistics. We also cal¬ 
culated an FPF of 2.7xl0” n (SNR~7) and 1.3xl0" 4 
(SNR~4) for b and c in F127M filter assuming a small 
number of resolution elements usin g the Student’s t-test 
relation presented in Mawet et al. (2014). 


4. RESULTS 

Contrast curves for the data in the F127M filter are 
shown in Fig. [2j the figure indicates the contrast achieved 
using different analysis techniques and the contrast in 
the global-KLIP and small-zone KLIP reductions. The 
global-KLIP reduction is not used in our analysis, we 
present this contrast curve estimate to guide future pos¬ 
sible studies at high-contrast with HST/WFC3-IR. The 
measured VEGAMAG photometry for the detections and la- 
limits is presented in Table [l] using revised zero-pointej/J 
The photometry of the planets is compared with brown 
dwarfs measured in the same HST filters in the color- 
magnitude diagram in Fig. [3] The positions of the plan¬ 
ets in Fig. [3j verifies that the flux measurement for the 
planets do not deviate unrealistically from the brown 
dwarfs. The spectral type derived for each planet is 
dependent on the contrast between J and water band 
fluxes being similar to that measured for brown dwarfs. 
Using a relation con yerting HST photometri c colors to 
a spectral type from Aberasturi et al. (|2014)f T °| we cal¬ 
culate a spectral type ot L9.5 ± U.5 for HR8799b and 
>L7 for HR8799c. The errors are the quadrature sum 
of the intrinsic scatter in the relation and the error in 
the photometric color. The larger role of clouds in gi¬ 
ant planets with temperatures comparable to cloud-free 
brown dwarfs complicates the conversion of the color- 
magnitude position into spectral types. The following 
two subsections model the combined HST photometry 

9 Revised ZP: F098M = 24.2209 mag, F127M = 23.7503 mag, 
F139M = 23.7679 mag 

10 SpT=1.56-6.25*(F127M-F139M) 
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with measurements from the literature at different wave¬ 
lengths to determine the atmospheric properties. 


TABLE 1 

HR8799BC ABSOLUTE PHOTOMETRY 


4.1. HR8799b 

The HST data in F139M provide the first detection 
of the planet b in the water absorption region at 1.4/mr. 
Additionally, the photometric points in F098M and the 
peak of the J band (F127M) are consistent with mea¬ 
surements taken ac ross similar wavelengths with ground 
based instruments (Currie c t al.|20lT Marois et al.|2008' 


Oppenheimer et al. 20fSb We combine the new HST 
data with photometric data available in the literature 
and fit the fluxes to sy nthetic model spectra from Bar¬ 
man et al. (2011 2015) (cited as B15 in the rest of the 
pap"ef). The models include clouds located at the inter¬ 
section of the pressure-temperature profile and the chem¬ 
ical equilibrium condensation curve with a parameterized 
thickness but a homogeneous distribution in latitude and 
longitude across the planet. The models also include non¬ 
equilibrium chemistry for all important C, N and O bear¬ 
ing molecules, as well as updated line lists for CH4 and 
NH3 (see B15 for details). The particle size-distribution 
is centered on 5/zm following a log-normal distribution. 

Figured] (top) shows two fits for the HR8799b photom¬ 
etry, the oest fitting model to the F13 9M ( black line) 


and the best fit to the (4.05] point (Currie et al. 2014 
green line). No individual model was able to lit the pair 
of neighboring fluxes at either F127M and F139M or L' 
and [4.05]. The discrepancy may indicate missing phys¬ 
ical processes in the model and the data is complicated 
by the non-contemporaneous measurements. Although 
the [4.05] band covers the Br-ct line, accretion is unlikely 
due to t he lack of Br-y emission in the K-band spectrum 
( Barman et al.|20TT I. The ratio of F139M flux to F127M 
flux may be useful for estimating the combined effects 
of cloud size distribution, cloud thickness and coverage. 
However, effective temperature, gravity and composition 
will need to be better constrained to avoid degeneracies 
in the model fitting and better agreement with evolution- 
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Fig. 3.— HST near-infrared color-magnitude diagram. The 
HR8799bc planets (red star, shaded region) are plotted with field 
L (blue circle s! and T (green circles') brown dwarfs from recent 
HST studies (|Aberasturi et al.|2014| |Apai et al.|2013| |Yang et al. | 
2014]. The black points are background held targets! The red 
rectangle for c represents our upper-limit estimate. 


Planet 

HST Photometry 

F098M F127M a F139M 

HR8799b 

HR8799c 

16.90 ± 0.18 16.20 ± 0.12 17.36 ± 0.26 

>15.72 15.38 ± 0.17 >16.16 


a P1640 spectrum through the F127M filter gives, 
HR8799b = 16.16i0.27 and HR8799c = 15.17i0.08. 


ary models. 

Fitting a combination of near-IR spectroscopy and IR 
photometry, B15 found a best matching model with T e fj 
= 1000K and log(g) = 3.5, with potentially subsolar 
water abundance and enhanced C/O, however, the un¬ 
certain surface gravity results in a wide spread of val¬ 
ues including solar metallicity. In this paper, the cooler 
best-fit model has solar abundances and we find T e ff = 
1000-1100 K and log(g) = 3.0-4.5, indicating consistency 
between photometry-based fits and those including spec¬ 
troscopy. Our measured water band flux, relative to J- 
band, however, is brighter than predicted by all of the 
non-solar models explored by B15, suggesting that per¬ 
haps clouds play an important role. 


4.2. HR8799c 

For HR8799c we measured the photometry in the 
F127M band and upper limits were determined for the 
F098M and F139M filters. The F127M flux is consistent 
with values obtained by integrating the s pect rum from 
PI640 over the same wavelength range (Oppenheimer] 
et al. 2013), but su bstantially lower th an the reported 
J-band photometry ( Marois et al.|[2008 ). The new HST 
data were combined with similar and longer wavelength 
fluxes from the literature to construct the spectral energy 
distribution given in Fig. [4] (bottom). Plotted with the 
data is a model with T ef j = HOOK, log(g) = 3.5 and ele¬ 


vated C/O t hat fit the K-band spectrum from Konopacky 
et al. (2013). The fits agree well with the HST photome- 


try and the longer wavelength data with the exception of 
a significant difference with the earliest measurement of 
the J-band photometry and an ^2.5cr difference with the 
[4.05] data point. Our models for HR8799c with homoge¬ 
neous cloud model agree, within the uncertainties, with 
the HST photometry and longer wavelength data out to 
~3.5^im with the exception of the first epoch of J-band 
photometrjrJ The difficulty in fitting all available pho¬ 
tometry witli reported uncertainties persists even when 
using patchy cloud models formed by combining clouds 
at different tempera tures (|Skemer et al.|2014|) or opacity 
(Currie et al. 2014) and the HST F1271V1 measurement 


is more consistent with the patchy cloud models. Our 
attempts with a linear combination of models did not 
provide a better fit than the homogeneous cloud model 
with higher C/O ratio presented in Fig. [4] (bottom). 


5. DISCUSSION 

Among the HST filters, the F127M is the most anal¬ 
ogous to a ground-based filter and can be directly com¬ 
pared to previous results reported in the J-band. The 
HR8799b F127M photometry is consistent with previous 


11 Improved analysis techniques and stellar variability cannot 
account for the factor of two difference in flux 
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Fig. 4.— Top: Model fits to HR8799b data. In the figure the red circles are HST photometry, and the blue squares are ground-based 



Integrating the HR8799c flux calibrated P1640 spectrum 
through the F127M filter matches the HST photometry 
to within 2 -a. A potential solution for this discrepancy 
might be intrinsic photometric variability caused by het¬ 
erogeneous cloud layers, which we find is not required 
in our model fits. The early groundbased photometry 
might also have suffered from calibration issues which 
combined by the intrinsic variability of the star might 
explain some of the difference in flux measurement. 

Efforts to match synthetic spectra to the ensemble of 
photometric data for HR8799b result in high effective 
temperatures and corresponding radii that are smaller 
by ~50% (0.69 - 0.92 R,i up ), than predicted by theoret¬ 
ical brown dwarf and giant planet cooling tracks (see 
Marley et al.| |2012 for a summary). Difficulty find¬ 


ing a model spectrum that simultaneously matches the 
near and mid-IR photometry could be due to the non- 
contemporaneous nature of the observations or perhaps 
the model-observation inconsistencies at multiple bands 


demonstrates that the atmospheres of both b and c are 
not enshrouded in high altitude hazes or clouds many 
pressure scale hei ghts thick, important for many tran¬ 
siting exoplanets (Kreidberg et al. 2014). Nonetheless, 
clouds are important in shaping the overall SED of the 
planets. Both best-matching models plotted in Fig. [4] 
(top panel) have clouds composed primarily of Iron and 
Magnesium-Silcate grains, located in the near-infrared 
photosphere. The cooler model has a cloud located at 
Pgas ~ 1 bar and extending upward 1 pressure scale- 
height. The warmer (and higher gravity) model has a 
cloud base near 10 bar, extending upward 2 pressure scale 
heights. Using single models to reproduce the observa¬ 
tions assumes global cloud coverage which is probably an 
overestimation. 

With spectral energy distributions well sampled obser- 
vationally and with model spectra that match reasonably 
well, the bolometric luminosities of both planets can be 
estimated, and we determine Lboi values of -5.1±0.1 for 
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b and -4.7±0.1 for c consistent with Marois et al. (2008). 


6. CONCLUSIONS 

In this study, we have demonstrated that WFC3 is 
capable of investigating the atmospheres of planets re¬ 
quiring high contrast at wavelengths either inaccessible 
from the ground or at which the extreme adaptive optics 
systems perform poorly. With its photometric stability, 
HST provides a valuable resource to explore exoplanet 
atmospheres and will enable measurements such as the 
search for atmospheric variability. The data from the 
current program can be used to build a PSF library for 
future high contrast imaging programs. Based on ex¬ 
perience from the current program future studies would 
benefit from obtaining more exposures rather than per¬ 
forming the time-consuming in-orbit roll. 

In this study, we were able to detect HR8799b in three 
WFC3 filters. In this work and in other studies, match¬ 
ing all of the data simultaneously has not been possi¬ 


ble. In particular, some of the shortest wavelength data 
favor warmer conditions while some of the longer wave¬ 
length IR data favor cooler conditions, possibly caused 
by clouds. For HR8799c, we find very good agreement 
between the model fit with the F127M, resolving the long 
standing difficulty of fitting the previously reported J- 
band photometry. A possible explanation for the F127M 
and J-band flux difference is variability. Future spec- 
trophotometric monitoring of HR8799b and HR8799c in 
the near and mid-IR might resolve whether the planets 
are variable over the duration of a single rotation period. 


Based on observations made with the Hubble Space 
Telescope, associated with program #12511. This re¬ 
search was supported in part by NASA cooperative 
agreements HST-GO-12511.04, HST-AR-12652.01, NSF 
AST-1411868. 
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